Charge and spin addition energies of one dimensional quantum dot. 
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We derive the effective action for a one dimensional electron island formed between a double 
barrier in a single channel quantum wire including the electron spin. Current and energy addition 
terms corresponding to charge and spin are identified. The influence of the range and the strength 
of the electron interaction and other system parameters on the charge and spin addition energies, 
and on the excitation spectra of the modes confined within the island is studied. We find by 
comparison with experiment that spin excitations in addition to non-zero range of the interaction 
and inhomogeneity effects are important for understanding the electron transport through one 
dimensional quantum islands in cleaved-edge-overgrowth systems. 

PACS numbers: 7f.10.Pm, 73.23.-b, 23.23.Hk 



I. INTRODUCTION 

Recently, considerable progress in the science and tech- 
nology of semiconductor nanostructures has been made. 
The experimental realization of one-dimensional (ID) 
quantum wires has opened new routes to investigating 
the influences of interactions and impurities on electron 
transport at low temperatures where quantum effects 
dominate In spite of numerous theoretical results 

obtained for the Luttinger liquid model , which have 
been used for interpreting the data, complete physical 
understanding of the experiments has not been achieved 
yet. It is therefore useful to study microscopically the 
combined effects of interactions and impurities on ID 
quantum transport within more realistic models and to 
compare the results with experimental data. 

The Luttinger liquid model allows to take into account 
interactions between Fermions in ID exactly. The linear 
conductance of a clean, infinitely long (spin degenerate) 
Luttinger liquid has been predicted to be G = 2e 2 g p /h, 
the universal conductance renormalized by an interac- 
tion constant g p 0. However, DC transport experi- 
ments on quasi- ID semiconductor quantum wires at rel- 
atively high temperatures showed only quantization in 
terms of the universal conductance G = 2e 2 /h. This 
has been interpreted either in terms of a compensat- 
ing, self-consistent rcnormalization of the external elec- 
tric driving field fsHlOfl or using non-interacting leads 
connected to the quantum wire [p]-[l3| which eventually 
determine the conductance. However, at lower temper- 
atures, temperature-dependent deviations from the uni- 
versal conductance steps have been found 
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When the electron density in the GaAs/AlGaAs-quan- 
tum wires fabricated by using the cleaved-edge-over- 
growth technique is decreased by applying a voltage to an 
external gate, eventually even the lowest electronic sub- 
band can be depopulated ITH] and the region of Coulomb 
blockade reached. Here, the mean electron density is suf- 
ficiently low, such that only very few maxima of the ran- 
dom potential of the impurities are higher than the Fermi 
level. A ID quantum island can be formed between two 
potential maxima in such a wire. Electron transport is 
then dominated by charging effects |i~6|-|lg|] . At tempera- 
tures lower than the charging energy, the linear conduc- 
tance shows discrete peaks corresponding to transferring 
exactly one electron through the quantum island. 

In this region, it has been detected that the tempera- 
ture dependence of the intrinsic width Y of several con- 
ductance peaks is modified by the correlations between 
the electrons and, instead of being independent of T, 
shows non-analytical power-law behaviour (A > 0) 
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(k < 0). This has been attributed to weak scattering at 
the random impurity potential in a Luttinger liquid Jl4| . 



Such a behaviour has been predicted in the sequential 
tunneling regime for resonant transport through a quan- 
tum dot coupled to a spinless Luttinger liquid |l9f| and 
for a double barrier in a Luttinger liquid p0pl| ]. 

In non-linear transport, collective excited states of the 
electrons in the island can contribute p2 23 1. The low- 
temperature current-voltage characteristics shows fine 
structure within the "Coulomb steps" which reflects the 
excitation spectrum of the electrons in the quantum dot. 
Especially in ID, strong effects of the spin of the corre- 
lated electrons have been predicted by using a sequen- 
tial tunneling model combined with the exact spectrum 
of eigenvalues and the corresponding states p4[ ]. In the 
above mentioned experiments |l5|, contributions of the 
excited states in the ID island have been found. In addi- 
tion, one of the corresponding peaks in the differential 
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non-linear conductance showed an interaction-induced 
non-analytical temperature dependence consistent with 
the above Luttinger liquid behaviour (||). Thus, correla- 
tions beyond the phenomenological charging model seem 
to be present in these semiconductor quantum wires. 
This is also indicated by recent results of Raman scat- 
tering experiments, which strongly suggest that semicon- 
ductor quantum wires are very probably dominated by 
non- Fermi liquid behaviour p5[ . In view of the existing 
broad theoretical understanding of ID non-Fermi liquids 
p6[ , it is highly desirable to provide quantitative results 
for charging and correlation phenomena in these systems 
by using reasonably realistic and controllable models, in- 
cluding non-zero range interaction between the electrons 
as well as the potential of the impurities. 

Efforts in this direction have been made by studying 
transport through a single barrier in a Luttinger liquid in- 
cluding long range interactions |^ |30| . Charging in the 
presence of two impurities has been studied for spinless 
electrons with interaction range much smaller than the 
distance between the two impurities ]3]J| . The crossover 
from a Luttinger liquid with long-range interaction to a 
Wigner crystal in the presence of two symmetrically ar- 
ranged potential barriers has been investigated |52| . 

In the present paper, we establish the fundamentals 
of a transport theory for ID quantum dots embedded in 
a ID electron system of correlated electrons with spin 
which turns out to be considerably different from the 
spinless case. Especially, we concentrate here on the en- 
ergies for charge and spin additions to the electron island 
formed between two strong impurities in a single-channel 
quantum wire of diameter d with an electron interaction 
of non-zero range. The quantum wire is described by 
using the Luttinger liquid model with spin ]26| . The 
bias electric field is assumed to have an arbitrary spa- 
tial shape. The interaction is assumed as a 3D Coulomb 
potential, which we project to ID by using the Gaus- 
sian wave function associated with the lowest state in a 
parabolic, cylindrically symmetric confinement potential. 
Screening is introduced by an infinite metallic plane at a 
distance D parallel to the wire. For comparison, we com- 
ment also on results obtained for a projected screened 3D 
Coulomb interaction. 

With the imaginary-time path-integral method, we de- 
termine the effective action in which charging and trans- 
port contributions are identified. The strength and the 
range of the interaction, Va and D, respectively, and the 
length of the ID quantum island - the distance between 
two delta-function-like potential barriers a - determine 
characteristic energy scales related to the charge and spin 
excitations that appear in the model. We find that the 
characteristic energy related to the addition of charge to 
the island - the charging energy - always increases with 
increasing interaction range until it saturates when the 
latter exceeds considerably the length of the island. On 
the other hand, the energy related to the addition of a 



spin - the spin addition energy - is approximately the 
same as that of non-interacting particles induced by the 
Pauli principle. This is due to the smallness of the (short- 
ranged) exchange interaction. 

We observe that the spatial shape of the electric field 
enters the effective action. The corresponding term can 
be interpreted as simulating the influence of the voltage 
at a gate electrode in experimental realizations. 

Using our results, we have analyzed recent experimen- 
tal data jl5| . We attempted to deduce the parameters of 
the quantum wire used in the experiment quantitatively. 
We found that it is not possible to identify a region of pa- 
rameters which is consistent with all of the experimental 
findings. We conclude that depending on the quantity 
considered, the frequency spectra of the collective exci- 
tations or the temperature dependence of the peaks in 
the differential conductance, experiments probe the in- 
teraction in different regions of the ID sample, namely 
within the island itself or within the whole quantum wire 
extending along the edge of the sample. This indicates 
that, together with spin and finite interaction range, in- 
homogeneity effects have to be taken into consideration 
for understanding the quantum transport in these ID 
semiconductor systems. 

The paper is organized as follows. Section || briefly 
describes the model. In Section III, the effective action 
is provided and discussed. In Section [lV| we comment 
on the conditions current transport. Quantitative results 
for the charging energy are given in Section |y|. Finally, 
we analyze experimental data in Section [V]]. 



II. THE MODEL 

A. The Hamiltonian 

We use the bosonization technique ||^] for inter- 
acting electrons with spin in ID. The quantum dot is 
described by a double barrier consisting of two delta- 
function potentials V, 5(x — Xi) at Xi(i = 1,2) where 
x\ < x 2 - An external electric field £{x,t) — —d x U(x,t) 
is assumed to induce transport. The Hamiltonian is 



H — Hq + Hb + Hh 



(3) 



The first term describes the interacting electrons 
within the Luttinger liquid model |f|,||, 



Hn = 



hv-f, 



dx 



n» + {d x v p {x)f 



dx l dx' d x ti p (x)V{x - x') d x S p (x') 
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dx 
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where v„ is defined below and 
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with the exchange interaction matrix element r] cx = 
V (2fcp) /2ttTivy . In Eq. (^), the electrons are represented 
by conjugate bosonic fields n p , -d p , and Ii ai d a associated 
with the collective charge and spin density excitations, 
respectively. The system length is L and the Fermi ve- 
locity v-p. The interaction, V(x — x'), is a projection of 
a modified 3D Coulomb interaction onto the x-direction 
(see below). It has the Fourier transform V(q) which is 
the dominant quantity in the dispersion relation of the 
charge excitations fl33|] 



u P {q) = VF \q\ 



Vex) 



2V(q) 



1/2 



v P {q)\g\ 



The dispersion relation of the spin excitations, 

1 1/2 _ , , 



Ucr(q) = v F \q\ [l 



Vex} 



(6) 



(7) 



contains as the dominant part the exchange interaction 
which is generally very small compared with V(q) at 
small q flj^S^. This implies that the group velocity of 
the spin excitations, v a is very close to the Fermi velocity, 
g a ~ 1. In the following, rj cx will therefore be omitted 
whenever it is only a small correction. 

The contribution of the two localized impurities is 



Hb = Pa 



i=1.2 



Vi cos [2kpXi 



x cos [v 27r^ CT (a;i)] 
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where po = 2k-p/n is the mean electron density. Equa- 
tion (g) is the potential energy of the impurities corre- 
sponding to the number density of charges 



p(x) = p ] {x)+p i {x) = p + \j-d x d p {x) 



in <M), besides the above Hb- As this represents a for- 
ward scattering process, it can be eliminated by a unitary 
transformation and will not be considered any further. 

On the other hand, if a potential is slowly varying on 
the scale fcp 1 its dominant contribution is mainly due to 
the long- wave length part of the density (^). We assume 
this to be the case for the bias electric field. The corre- 
sponding term in the Hamiltonian is, with the elementary 
charge e(> 0) 



H„ = 



-e-W- / dxU{x,t)d x d p {x). 



(11) 



The presence of the two localized impurity naturally 
separates the charge and spin degrees of freedom at 
"bulk" positions x ^ x\^x% from those at the barriers. 
It is useful to introduce symmetric and antisymmetric 
variables for particle (y = p) and spin densities (y = a) 
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The quantity N~ is associated with the fluctuations of 
the particle number within the island as compared to 
the mean particle number no = po(x2 — x±). The corre- 



sponding excess charge is Q 



-eN ; N a represents 



the z-component of the fluctuation of the number of 
spins within the island corresponding to a change of spin 
N~ /2. The numbers of imbalanced particles and spins 
between left and right leads are represented by N$ . The 
DC current-voltage characteristic can then be evaluated 
by considering the stationary limit of the charge transfer 
through the dot in the presence of an external voltage 



2 t 



(13) 



The brackets (• • •) include a thermal average over the 
collective excitations at x ^ x\,X2, and a statistical av- 
erage performed with the reduced density matrix for the 
degrees of freedom at x = x\, x 2 - 



+p cos [2k F x + v2tti9 p (x)] cos [y27n? CT (i)]. (9) 



B. The Interaction Energy 



The second term in Eq. (g) accounts for the slowly vary- 
ing part of the charge fluctuations. The third term 
represents the charge density wave involved in the 2fcp 
backscattering interference between left and right mov- 
ing electrons. It also couples the charge with the long 
wave length part of the density of spins, 



Pt( x )-Pl( x ) 



- d x d a (x) , 



(10) 



which is considered here with respect to a zero mean 
value. When calculating the energy of the impurities, 
one obtains also a contribution due to the second term 



We consider in the following two models which can be 
used for discussing experimental results quantitatively. 



1. Model 1. 

Because of gates and surrounding charges in the ex- 
periments done on quantum wire, screening is always 
present for the interaction between the electrons which 
normally interact via the 3D Coulomb interaction. Al- 
though the geometry of the gates and the wire are cer- 
tainly more difficult in the experiment |l5|, we assume 
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in the following that we can describe the screening as 
being solely due to an infinite metallic plane at a dis- 
tance D in the y-direction parallel to the [x, z)-plane. 
The quantum wire is assumed in the x-direction. The 
effective electron-electron interaction energy can then be 
calculated by using the method of image charges 



V(r-r') 



V 



y/(x - x') 2 + {z- z' f + (y + y'- 2D) 2 



(14) 



with Vo — e 2 /47reoe and the dielectric constant eye. 

We assume for simplicity a parabolic confinement of 
the electrons perpendicular to the wire. The effective 
interaction in the lowest subband can then be obtained 
from (fbi]) by projecting with a normalized Gaussian con- 
finement wave function corresponding to the diameter d 
of the wire. The Fourier transform of the resulting effec- 
tive interaction potential is for D 3> d 



V{q) = Vo 



3 w Ei W 



2K (2Dq) 



(15) 



where Ei(z) is the exponential integral and Ko(z) the 
modified Bessel function This expression shows 

explicitely how the gate screens the Coulomb interac- 
tion. In the limit qd — > 0, one obtains the finite value 
V{q -> 0) = 2V Q [7/2 + log(2£)/d)] (7 = 0.57722 Euler 
constant). This implies a finite interaction constant 



9p 



i'F 



v P (q -> 0) 



1 + rpf + 2?/ log 



2D 
~d 



-1/2 



(16) 



where r\ = ZVo/ithvp- The ID equivalent of the un- 
screened Coulomb interaction is obtained for D — ► 00. 



2. Model 2. 

We also consider the screened Coulomb interaction 



V(f-r') = Vo- 



ir — r 



(17) 



with a phenomenological screening length a -1 . The 
Fourier transform of its Gaussian projection to ID is |l(]] 



V(q) = Voe^/^+^E, {^-[q 2 + a 2 
and the corresponding interaction constant 



9p = 



1 - 777 + 2?? log 



-1/2 



(18) 



(19) 



III. THE EFFECTIVE ACTION 

In order to evaluate the current-voltage characteris- 
tic one has to perform a thermal average over the "bulk 
modes" at x 7^ Xi,X2- This can be done with the 
imaginary-time path integral method pjJ. The result 
of the integration is an effective action S e s which de- 
pends only on the four variables defined in (|l^). In the 
continuous limit [L — > 00), with the inverse temperature 
(3 = 1/kftT one obtains 



S eS {Nf,Nj 

+ E E 

r=± u—p,(T 



hp 



dr H B [N±,N±] 



K/3 rh/3 



drdr' NI{t)KI{t-t')NI{t') 



H0 

.„.„. drN r p (r)C r (r) 
/o 



(20) 



The Fourier transforms, at Matsubara frequencies lu h — 
2irn/Ti(3 1 of the dissipative kernels -ftTjf (r) and of the ef- 
fective "forces" £ ± (t) are determined by the dispersion 
relations of the collective modes (||) and (^) , respectively, 



1 ± cos[g(xi - x 2 )} 
^n+uKq) 



. (21) 



' ' ' — ~j~T K p k ( LJn ) I dx£{x,u n ) 



dq 



cos[q(x — X2)] ± cos[q(x — xi)] 



(22) 



Both, and £ ± contain the collective bulk modes 
which introduce the interaction effects to be described 
below. First of all, we note that K+(uj n -> 0) = {HUH. 
On the other hand, K ~ (u> n — > 0) ^ 0. The latter describe 
the costs in energy for changing the numbers of charges 
and/or spins on the island between the potential barriers. 
The corresponding Euclidean action is 



S [N~,N- 



E 

U—p,(T 



2 



h/3 



dr (N~ 



with the characteristic energies 



E v = 2K-{w n ^Q) {v = p,a). 



(23) 



(24) 



For v — p, this corresponds to the charging energy that is 
supplied/gained, in order to transfer/remove one charge 
to/from the island as compared with the mean value, 
Np = ±1. Correspondingly, for v = a, the spin addition 
energy E a is needed/gained in order to change the spin by 
exactly ±1/2. The Coulomb interaction in the dispersion 
relation of the charge excitations, increases considerably 
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the charging energy E p , in comparison with the spin ad- 
dition energy which is only influenced by the (small) 
exchange interaction. We always expect E p > E a . 

The frequency dependent parts of the kernels describe 
the dynamical effects of the external leads and of the 
correlated excited states in the dot. Their influence is 
described by spectral densities J^(cj) related via analytic 
continuation to the imaginary-time kernels 20 3q] 



J±H = — hnK±(u r . 



(25) 



Due to the non-zero range of the interaction, analytic 
expressions for these densities are not available. However, 
one can always extract their limits for lo — > 0, 



J±(w^0) 



±9v 



(26) 



where A~ = g 4 p (E p /E ) 2 (E = Knv F /2a), 
three other combinations of indices A^; = 1 



and for the 
This limit 

describes the dissipative influence of the low-frequency 
charge and spin excitations in the external leads, x < x\ 
and x > x%. It holds also for finite frequencies. How- 
ever, these must be smaller than the frequency scale cor- 
responding to the range of the interaction, and smaller 
than the characteristic excitation energy of the correlated 
electrons in the dot. 

Let us now discuss the driving forces. In general, 
C^{t) depend in a quite complicated way on the dis- 
persion of the collective modes and on the shape of the 
electric field. We focus on the DC limit where it is suf- 
ficient to evaluate the Fourier components for uj n — ► 0. 
In this case, the quantity £ + (t), which acts on the total 
transmitted charge, depends only on the integral of the 
time independent electric field over the entire system, the 
source-drain voltage U = J_ dx£ (x), 



eU 



(27) 



Since C + is the part of the effective force that generates 
the current transport, this result generalizes the one ob- 
tained previously for only one impurity J35|. It can be 
easily derived from Eq. (E2h by using the relation 



e 2 VF 
Im 2 



dq 



u> n (l ± cos qx) 
= a (0,u) n ) ± a (x,uj n ) 



(28) 



Here, <7o(x, uj n ) is the frequency dependent non-local con- 
ductivity of the Luttinger liquid per spin channel [j35f , 
with the DC limit o-q(x,0) = g p e 2 /h. 

On the other hand, £~{t) acts on the excess charge on 
the island, it does not generate a current. It depends on 
the spatial shape of the electric field and can formally be 



written in terms of the total charge Qs accumulated be- 
tween the points x\ and xi in the absence of the barriers 
as a consequence of the presence of the DC electric field 



C~(r) = 



EpQe 



(29) 



where the charge is given by 
Q £ =2 / dx'S(x') 



x lim 



a (xi - x' , -iuj) - o-q(x2 - x 1 , -ioj) 



(30) 



Equivalently, this can also be understood in terms of 
addition energies. By introducing explicitly in Eqs. ( p^ , 
pi| ) the dependence on the interval considered when eval- 
uating the addition energies, one finds 



c ( r ) = 2 E p( Xx ~ X2 > 



dx£{x) 



E p {x-xx) E p {x-x 2 ) 



(31) 



It is reasonable to assume £2,1 = ±a/2. If the effec- 
tive electric field has inversion symmetry, Cr vanishes. 
Without inversion symmetry, the electric field generates 
an effective charge on the island which will influence the 
total current via coupling between and N~ due to 
the impurity Hamiltonian Hb ■ Physically, this externally 
induced charge may be thought of as being generated by 
a voltage Vq applied to an external gate which electro- 
statically influences the charge on the island. Thus, the 
above Eq. ( p9| ) can be interpreted as a term represent- 
ing the effect of the gate voltage in the phenomenological 
theory of Coulomb blockade. 



IV. CONDITIONS FOR CURRENT TRANSPORT 

So far, we have discussed the influence of the spin and 
charge bulk modes on the four variables JV^ that describe 
the quantum dot. In order to calculate explicitly the 
electrical current one has to solve the equations of motion 
for the N„ . This is beyond the scope of the present work 
and will be discussed elsewhere. Here, we briefly describe 
the results that are needed in the following. 

For barriers much higher then the charging energy, 
Vi E p , the dynamics is dominated by tunnel- 

ing events connecting the minima of Hb in the 4D 
(7V+,7V / 7,iV+,7V ( 7)-space 0. For equal barriers, V x = 
V2 = V, the impurity Hamiltonian is 



H B = PoV 



nN+ 
cos — - 1 — COS 



ttN+ 7r(n 
— - — cos 



TV; 



■ cos ■ 



■ sin ■ 



. ir(n + N-) . nN a 
■ sin — sin — — 



(32) 
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The transitions between the minima of this function of 
four variables correspond to different physical processes 
of transferring electrons from one side to the other of the 
quantum dot. At very low temperature, the dominant 
processes transfer the electron coherently through the 
dot. In particular, when the number of particles in the is- 
land is an odd integer there is spin degeneracy, N~ = ±1. 
The island acts as a localized magnetic impurity, similar 
as in the Kondo effect j7|. This leads to a Kondo reso- 
nance in the transport through the island fl37p8] ]. 

On the other hand, if the temperature is higher than 
the tunneling rate through the single barrier, the dom- 
inant processes are sequential tunneling events |]l9| , pT| . 
The transfer of charge occurs via uncorrelated single- 
electron hops into, and out of the island, associated with 
corresponding changes in the total spin. This is precisely 
the regime that recently has become accessible by using 
cleaved-edge-overgrowth quantum wires . In this re- 
gion, one has to consider minima which correspond to 
pairs of even or odd N~ and N~ . The dominant trans- 
port processes are those which connect minima via tran- 
sitions N~ — > N~ ± 1 associated with changes of the 
spin N~ — » N~ ± 1. For each of these processes also the 
external charge and the spin change by iV+ — > N+ ± 1. 

The degeneracy of these minima is lifted by the charge 
and spin addition energies E p and E a which force the 
system to select favourable charge and spin states in the 
island. These selections become essential at low tem- 
peratures, k^T < E Vl when current can flow only under 
resonant conditions. The latter can be achieved in experi- 
ment by tuning external parameters, like the source-drain 
voltage or the gate voltage, in order to create degenerate 
charge states in the island. 



A. Linear Transport 

In the linear regime, U — > for T = 0, starting with 
the island occupied by n electrons, we expect that an- 
other electron can enter and leave only if the difference 
between the ground state energies of n + 1 and n elec- 
trons is aligned with the chemical potential of the ex- 
ternal semi- infinite Luttinger systems. The ground state 
of an even number of electrons in the ID island has the 
total spin 0. On the other hand, the ground state of an 
odd number of electrons can be assumed to have the spin 
N~ = ±1 p4|,p9|. This implies the resonance condition 



U (n + 1, ±s„+i) ~U(n, ±s n ) = 



(33) 



with U (n, ±s n ) the ground state energies with n particles 
and total spins s„ = (n even) or s„ = 1/2 (n odd), 
respectively. 

With the above charge, spin and external gate terms, 
Eqs. (E4) and (E9), these conditions become 



E p (n-n -n G + ^j+(-l) n ^ = 0. (34) 

The variable uq = eVQ5/E p represents the number of 
induced particles due to the coupling to a gate at which 
the voltage Vg is applied, with a proportionality factor 
S which can be determined experimentally. The zero of 
energy has been assumed to be given by the external 
chemical potential in Eq. (|3^). 

From the above expression one can see that the dis- 
tance of the peaks of the linear conductance when 
changing the gate voltage are given by AVq = (E p + 
(—l) n E a )/e5 w E p /e5 since E a <C E p . Having indepen- 
dent information on 8 one can extract the value of the 
charging energy E p (in principle also for the spin addition 
energy E a ) from the experimental data. 

In order to evaluate the current as a function of tem- 
perature and/or source-drain voltage, one needs to con- 
sider the behavior of the spectral densities given in 
Eq. (|2^). In the sequential tunneling regime, one can 
show that only their sum enters the transport Elt] 



■%) = E E J »- 



(35) 



The frequency behavior of this determines the current- 
voltage characteristics both in the linear and in the non- 
linear regimes (see below). In particular, it determines 
the exponent of the power-law dependencies of the cur- 
rent as a function of temperature and/or the bias voltage. 

For temperatures lower than the excitation energy in 
the quantum dot, the low-frequency behavior of the spec- 
tral density is mainly determined by the charge and spin 
excitations in the leads. Thus, the power-law dependence 
of the conductance peaks is dominated by the interaction 
in the regions of the quantum wire outside of the electron 
island. From Eqs. (p5|,p6|) one obtains 



J(w) W JlcadsM 



2 



1 



i + A: 



(36) 



Eq. ([36]) generalizes the results obtained previously 
for spinless electrons and zero-range interaction [ ]l9| , pl"| |, 
where J(cj) « u/g p . We conclude that the presence of 
the spin in the leads introduces an effective interaction 
strength 



1 

geff 



1 + A: 

f_ 

9p 



2 

9<y 



(37) 



which determines the exponents of the power laws. 

For example, it has been shown for spinless electrons 
p9| , pl| that the intrinsic width T(T) of the linear con- 
ductance Coulomb peak depends on the temperature. 
For low temperatures this has been found to be given by 
T(T) cx T 1 / 9 "^ 1 . Correspondingly, with spin, one finds 
instead 



r(T) cx T 



l/Se: 



(38) 
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B. Non-Linear Transport 



V. THE CHARGING ENERGY 



One can investigate the non-linear current-voltage 
characteristic by increasing the source drain voltage U. 
In this case, the current-voltage characteristic shows the 
Coulomb staircase associated with transitions between 
successive ground states of the electrons in the quantum 
dot. In addition, fine structure appears which reflects the 
excitation spectra of the correlated electrons. They can 
be due to collective charge and spin modes or to partic- 
ular polarizations of the spins, N~ =/= 0, ±1. The former, 
for a voltage drop U at the potential barriers (fixing the 
chemical potential in one of the leads to be equal to that 
in the dot), has a maximum periodicity 



U = fx (n + 1) - no{n) = E P + (-1) 



n+l 



E„ 



(39) 



where no(n) is the electrochemical potential of n elec- 
trons on the quantum island. 

In the present model, the possible particle-hole excita- 
tions are collective spin and charge modes confined within 
the island. In a completely isolated island, the corre- 
sponding energy spectrum would be discrete, uj p {q m ) and 
w ff (q m ), due to the discretization of the wave number 
Qm = trm/a associated with the confinement. Due to the 
coupling to the electrons in the leads via the interaction, 
these levels are broadened. In the following argument, 
we assume that this broadening is negligible. 

The screened Coulomb interaction causes a non-linear 
dispersion relation for the charge modes in the infinite 
Luttinger system. This leads to non-equidistant charge 
excitation energies in the quantum dot, 



Ae p (q m ) = H [oj p (q m+ i) - oj p (q m )] . 



(40) 



Their explicit values depend on the ratio between the 
distance a and the range of the interaction, D or a -1 . 

For values of a much larger than this range, the first 
excited charge modes are equidistant. They are given by 
the charge-sound velocity for q — > 0, v p — vp/g p , 



hlTVn 



hnvp 



2Eo 
9p 



(41) 



In the opposite limit, the non linear dispersion is already 
affecting strongly the first excitation. This results in a 
value smaller than Eq. (^|), for model 1, 



hirvp 



1 — 777 + 2?y log 



2a 
nd 



2Eo 



(42) 



For the spin excitations, the dispersion in the infinite 
system is linear and the discrete spectrum is equidistant 



2E - 



a vp 
with the spin mode velocity tv(~ vp). 



(43) 



The above discussion emphasizes the importance of 
charge and spin addition energies, and of the charge 
and spin excitation energies for the linear and non-linear 
transport properties. It is therefore crucial to evaluate 
these quantities microscopically and determine their de- 
pendencies on the model parameters, especially in view 
of comparisons with the experimental data of pi . 

We analyze now in more detail the charging and spin 
energies of the Luttinger island as a function of the pa- 
rameters of above Model 1. In addition to the interaction 
strength Vq, which contains as an essential quantity the 
dielectric constant e, we have the Fermi energy Ep and 
three geometrical parameters: the distance D between 
the ID system and the infinite metallic plane, the diam- 
eter d of the wire and the length a of the island. 

The spin addition energy can be easily evaluated from 
the simple dispersion relation Eq. (0) 



E a = . 

2a g a 



(44) 



It is the same as the addition energy of non-interacting 
particles which is related to the Pauli principle and due 
to the quasi-discretization of the spin energies in the dot. 

The charging energy E p is evaluated numerically from 
Eq. ( p4f ) as a function of a/d. It is shown in Fig. |l| for 
different ratios D j d and different interaction parameters 
Vo/d with d = lOnm and d = 20nm. Different Vo/d 
correspond to different e. Changing Ep between 2meV 
and 4 meV, values that cover the experimentally relevant 
region, changes the curves in Fig. |l| only less than 10%. 
For small a/d, the charging energy diverges. We consider 
only the region a/d > 1. Asymptotically, 



E ' 

El, 

E ' 



1 + 777 + 2i] log 



1 - 777 + 277 log ( — 



2D 



2 a 



1 



1 

72 



(D<o) 



(45) 



(D»a). (46) 



E p is only weakly influenced by D/d: changing D/d by a 
two orders of magnitude changes E p only by about 30%. 

The charging energy always increases with increas- 
ing interaction range D, because the cost of energy for 
putting additional electrons into the island increases. In- 
teraction ranges much larger than a do not change the 
charging energy because only the short-range part of the 
repulsion between the electrons contributes. Therefore, 
E p approaches the asymptotic value of Eq. ([hj|) (Fig. [l], 
curves for largest D/d). For strong Coulomb interaction 
(77 = 2Vq/itTivp ^> 1) and a^> d, this is 



E n 



27reena 



log 



e 
C 



(47) 
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Here, the C is the classical self-capacitance of a cylinder 
of the length a. The stronger the interaction (increasing 
Vo) the larger is E p , very similar to the behaviour of a 
classical capacitor with a dielectric, C oc s. 

Very similar results, with d/D replaced by ad, are ob- 
tained by using Model 2. Basically, this tells that screen- 
ing effects may well be described by a global parameter, 
irrespective of the underlying microscopic model. 

As mentioned earlier, the non-linear transport shows 
Coulomb steps with widths given approximately by E p , 
and fine structure due to the excited states. For a first 
estimate of the strength of the interaction is very useful to 
estimate how many excited states are present within the 
energy window given by E p . Figure shows the ratios 
between the energy E p and the first collective excited 
(spin or charge) state in the electron island Ae„ [y = a, p) 
for Model 1 as a function of the interaction constant g p . 
The change in g p is obtained changing the distance D 
from the gate, g p — > corresponds to D — - > oo. The ratio 
E p /Ae p is always smaller than E p /Ae a . 

We expect to observe within the window E p much more 
levels due to spin excitations than due to charge excita- 
tions; E p / Ae v always increases when g p decreases (in- 
creasing D). It saturates for D 3> a, 



3l 

Ae n 



1 



IlL 
Ae a 



In the opposite limit, D <C a, the behavior is described 
by the asymptotic expressions 



IlL 
Ae 



1 

2g~ P 



3l 

Ae 



Vp 1 



(49) 



VI. COMPARISON WITH EXPERIMENT 

In this section, we compare our results with the ex- 
perimental data obtained in p5[ . In that work, results 
of the temperature dependence of the intrinsic width of 
the conductance peaks in the region of Coulomb block- 
ade on quantum wires fabricated with the cleaved-edge- 
overgrowth technique have been reported. The data have 
been found to be consistent with the power law 



T(T) cx T 1/9 " 



(50) 



with g* ft 0.82 and g* ft 0.74 for a peak closer to the 
onset of the conductance and the next lower peak, re- 
spectively. In addition, information about the excited 
energy levels of correlated electrons have been obtained 
by measuring the non-linear current voltage characteris- 
tics. A minimum of five excited levels have been observed 
for a given electron number. Presumably, the number of 
excited levels is even larger since the majority of the ex- 
cited island states is generally not clearly observed in 
non-linear transport due to matrix element effects ECf]. 



The data have been analyzed previously by assuming 
that within the quantum wire, a quantum dot has been 
accidentally formed between two maxima of the random 
potential of impurities. The electron spin has been ne- 
glected. The parameters estimated from the experimen- 
tal setup are: length of wire L ft 5 /im; length of the 
electron island a ft 100 — 200 nm; (non-spherical) diam- 
eter of the wire d ft 10 — 25 nm; distance to the gate 
D ft 0.5 /im. The charging energy, as estimated from the 
distance between the conductance peaks, has been given 
as Eq ft 2.2 meV. Since the observed conductance peaks 
are equidistant within 10% we deduce that E„ <C E p . A 
rough estimate of the Fermi energy is Ep ft 3 meV. With 
these parameters, the value of the interaction constant 
has been estimated as g p ft 0.4 [|||, clearly inconsistent 
with the above g* determined from the temperature de- 
pendence of the conductance peaks. 

Using our above microscopic approach which includes 
the influence of the electron spin and finite range of 
the interaction, we confirm that there is a discrepancy, 
though the estimate for g p from the temperature depen- 
dence of the conductance peaks seems to be slightly im- 
proved. 

First, we estimate the length of the island a from the 
charging energy in Fig. [l]. As E p is relatively insensi- 
tive with respect to changes of D/d and Ep in the ex- 
perimentally relevant region of parameters, we assume 
D/d = 50 with E F = 3meV and d = 20nm. Further- 
more, Vo/d — 6.02 meV for e = 12. With the above 
2.2 meV we find a/d ft 16. This is consistent with the 
value given in p| within the uncertainties. 

By taking into account the spin and using an interac- 
tion of a finite range, g* has to be identified with the 
effective interaction g e f[ in Eq. (|37]). Since for the spin 
excitations g a ft 1, we find for g* = 0.82 and 0.74 a in- 
teraction constant g* ft 0.6 ±20% by solving Eq. (|3^ ) for 
E p (gp)/Eo and comparing with E p (g p ) /2E$ = E p /Ae a 
in Fig. | This would imply that D ft d (Eq. ©). On 
the other hand, by using Eq. (^) with the above param- 
eters, especially D ft 50d, we find a interaction constant 
g p ft 0.25 ± 20% depending only weakly on D/d. This 
is clearly not a quantitatively consistent result. The dis- 
crepancy can be reduced by changing e. Also, taking into 
account the exchange correction r] cx (g a ^ 1) gives some 
improvement. However, within reasonable values of all 
of the parameters, one is in any case forced to conclude 
9* P « 2 5p . 

Furthermore, with g* = 0.6, one reads from Fig. g 



E p /Ae 



Ep/Ae 



1. Thus, about 1 to 2 excited 



states corresponding to a given electron number should 
be observed in non-linear transport spectroscopy. This 
is also not consistent with the experiment. In order to 
observe a larger number of excited states, the interaction 
constant must be considerably smaller, g p < 0.3, such 
that the number of spin excited states within the given 
window of E p is increased (cf. Fig. 0). Taking into 
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account the exchange interaction [v a < vp) would ad- 
ditionally decrease Aeo- (and E a ) increase E p , and thus 
increase the number of excited states. 

Despite the above inconsistency, the energetically low- 
est of the excited states seen in the non-linear transport 
are predicted by our approach to correspond to spin exci- 
tations. Without the latter there is not at all any possi- 
bility to understand the results from the non-linear trans- 
port even when using parameters that give consistent re- 
sults for the interaction constants. Given this scenario, 
only the energetically highest experimentally observed 
state would correspond to a charge excitation. The latter 
could be identified, for instance, with the state denoted 
in (l^] as "strongly coupled excited state" . 

One might conclude that the part of above quantitative 
discrepancy is due to the Luttinger model being question- 
able for quantum wires in the extreme low-density limit. 
However, one can also conclude that the data obtained 
from non-linear transport spectroscopy have to be inter- 
preted by using a different interaction constant than that 
obtained from the temperature dependence of the trans- 
port. This is also supported by the theoretical derivation: 
the temperature behaviour of the conductance peaks is 
dominated by the excitations in the whole quantum wire, 
while the quasi-discrete excitation spectra are showing 
the interaction in the region of the quantum island. We 
are thus forced to conclude that inhomogeneity effects 
and the influence of the contacts are an important is- 
sue for the understanding of the correlations in the elec- 
tron transport in these cleaved-edge-overgrowth quantum 
wires, in addition to spin effects. 

VII. CONCLUSION 

In conclusion, we have derived the effective action for 
a quantum dot formed by a double barrier potential with 
a realistic long range interaction between the electrons, 
and, most importantly, including the electron spin. Wc 
have identified an effective driving force acting on the 
charge of the electrons, which turned out to be indepen- 
dent of the shape of the driving electric field. Mass terms 
originating from the dissipative degrees of freedom in the 
Luttinger liquid have been found for both the charge and 
the spin. They are identified with the charging energy in 
the case of the addition of charge and with a spin addi- 
tion energy for the spin. Also, an effective force has been 
found which, in the DC limit, is only non-zero when the 
driving electric field and/or the barriers are asymmetric. 
This influences the transport via the coupling of modes 
at the potential barriers. It can be attributed to charging 
of the electron island via an external gate. 

We have made an attempt to interpret quantitatively 
a recent experiment which includes linear as well as non- 
linear transport data in the region of Coulomb blockade. 
We find that even taking into account the electron spin 



it is impossible to consistently fit all of the experimental 
data with the same interaction parameter. It seems that 
different interaction strengths have to be used for ex- 
plaining the temperature dependence of the transport on 
the one hand, and the excitation spectrum of the quan- 
tum dot on the other. 

A possible conjecture to solve this puzzle is to assume 
that the temperature behaviour of the conductance peaks 
is dominated by the interaction among the electrons in 
the quantum wire along the edge of the whole sample, in 
contrast to the excitation spectra which feel the interac- 
tion strength near the quantum dot. This is qualitatively 
consistent with the physical origin of the temperature 
dependence of the conductance, namely the bulk modes 
outside the electron island. This conjecture would also 
explain why g* rs 0.6, as estimated from the temperature 
behaviour, is larger than that obtained from the observed 
number of the excited states in the island, g p < 0.3. The 
interaction is weaker in the leads due to screening in- 
duced by the contact regions with the attached 2D elec- 
tron gases while it can be assumed to be stronger within 
the island due to the absence of nearby 2D electron gases. 

The above interpretation also predicts that the en- 
ergetically lowest collective excitations in the quasi-lD 
quantum island observed in the experiment are very 
likely spin excitations. The energies of the lowest ob- 
served collective excited states are so small that it is hard 
to believe that they can be related to charge modes since 
v p > v a , given the above estimates for the range and 
strength of the interaction is of the correct order. This 
is consistent with a previous suggestion based on a semi- 
phenomenological model in which exact diagonalization 
and rate equations have been combined in order to de- 
scribe the non-linear transport through a ID quantum 
island containing only few electrons ]24| . 
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FIG. 1. Charging energy E p in meV for Model 1 as a 
function of a/d (Fermi energy Ef — 3meV, effective elec- 
tron mass m = 0.067mo). Top: d — lOnm; interaction 
strengths Vo/d = 24.08 meV (dashed-dotted), 12.04 meV (full 
lines,) 6.02 meV (dotted), corresponding to e = 6, 12, 24 
(equivalent to r\ — 2Vo/nhvF = 2.04, 1.02, 0.51, respec- 
tively), and D/d = 1000, 100, 10 (top to bottom). Bot- 
tom: d = 20nm; interaction strengths Vo/d = 12.04 meV 
(dashed-dotted), 6.02 meV (full lines,) 3.01 meV (dotted) and 
D/d = 500, 50, 5 (top to bottom). 
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FIG. 2. Ratios Ep/Ae^ (u — a, p; upper/lower curves for 
fixed a/d) as a function of the interaction constant g p for 
a/d = 20, 60, 500 with e = 12, d = 20 nm, m = 0.067m , 
9<y = 1- 
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